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investigated by using a method called as the generalized tortoise coordinate 
transformation. The location and the temperature of the event horizon of the 
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I. INTRODUCTION 

The fourth quarter of last century has witnessed various remarkable progress on several 
researches on black hole physics since Hawking's marked discovery • One of these aspects 
is to reveal the thermal properties of various black holes 0. Much efforts have been devoted 
to studying the thermal radiation of scalar fields or Dirac particles in the stationary axisym- 
metry black holes A popular method to determine the location and the temperature 

of the event horizon of an evaporating black hole is to calculate vacuum expectation value 
of the renormalized energy momentum tensor ||. But this method is very complicated, it 
gives only an approximate value of the location and that of the temperature. Thus it is of 
limited use and meets great difficulties in many cases. 

Zhao and Dai suggested a novel method of the generalized tortoise coordinate trans- 
formation (GTCT) which can give simultaneously the exact values both of the location and 
of the temperature of the event horizon of a non-stationary black hole. Basically, it is to 
reduce Klein-Gordon or Dirac equation in a black hole space-time to a standard wave equa- 
tion near the event horizon by means of generalizing the common tortoise-type coordinate 
r* = r + j- ln(r — r#) in a static or stationary space-time @,[|[5] (re is the surface gravity of 
the studied event horizon) to a similar form in a non-static or non-stationary space-time H 
by allowing that the location of the event horizon r# can be a function of the advanced time 
v = t + r* and/or of the angles 9, (p. A key point in introducing the tortoise coordinate r* is 
that one must base upon the symmetric property of the considered space-time. For exam- 
ples, the location of the event horizon is a constant (r# = 2M) in the Schwarzschild black 
hole, while it is a function of the advanced time (r# = rjj(v)) in a Vaidya-type space-time. 
The GTCT method has been applied to investigate thermal radiation of scalar particles in 
the case of non-stationary axisymmetric black holes 0. Hawking effect of Dirac particles in 
the non-static black holes has also been studied successfully in . 



However, it is very difficult to discuss on the evaporation of Dirac particles in a non- 
stationary axisymmetric black hole. The difficulty lies in the non-separability of the radial 



and angular variables for Chandrasekhar-Dirac equation [11 1 in a non- stationary axisymme- 



try space-time. In this paper, we try to tackle with the thermal radiation of Dirac particles 



in a non-stationary Kerr space-time ||12] , |13| . We consider the asymptotic behaviors of the 
first-order and second-order forms of Dirac equation near the event horizon. Using the re- 
lations between the first-order derivatives of Dirac spinorial components, we can eliminate 
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the crossing-terms of the first-order derivatives in the second-order equation and recast each 
second-order equation to a standard wave equation near the event horizon. The location 
and the temperature of the event horizon are just the same as that obtained in the case 
of the thermal radiation of Klein-Gordon scalar field in a non-stationary Kerr space-time, 
but the Fermionic spectrum of Dirac particles displays another new effect dependent on the 
interaction between the spin of Dirac particles and the angular momentum of black holes. 

Within the framework of Newman-Penrose formalism [14}|, we derive in Sec. II the explicit 
form of Dirac equation in the variable- mass Kerr space-time. In Sec. Ill, we introduce a 
GTCT according to the symmetry property of space-time and investigate the asymptotic 
behavior of the first-order Dirac equation near the event horizon. Then the equation that 
determines the location of the event horizon are inferred from the vanishing determinant of 
the coefficients of the first-order derivative terms. A crucial step of our treatment in Sec. IV 
is to use the relations between the first-order derivative terms to eliminate the crossing term 
of the first derivatives in the second-order Dirac equation near the event horizon. Then we 
adjust the temperature parameter k introduced in the GTCT so as to recast each second- 
order equation into a standard wave equation near the event horizon. In the meantime, we 
can get an exact expression of the Hawking temperature. Sec. V is devoted to separation of 
variables and to presenting the thermal radiation spectrum of Dirac particles according to 
the method of D amour- Ruffini-Sannan's Finally, we give a brief discussion about the 

spin- rotation effect. 



II. DIRAC EQUATION IN A NON-STATIONARY KERR SPACE-TIME 



The variable mass Kerr solution |T2|,|T3j can be written in the advanced Eddington- 
Finkelstein system as 

A - a 2 sin 2 9 



ds 2 = aV - Idvdr + 2a sin 2 9drdip - T,d9 z 

r 2 + a 2 -A . 2 (r 2 + a 2 ) 2 - Aa 2 sin 2 9 . 2 2 
+2 a sm 2 Odvdip - - sm 2 9 dip 1 , (1) 

where A = r 2 — 2M(v)r + a 2 , S = r 2 + a 2 cos 2 9 = p*p, p* = r + ia cos 9, p = r — ia cos 9, 
and v is standard advanced time. The mass M depends on the advanced time v, but the 
specific angular momentum a is a constant. 

We choose a complex null-tetrad {I, n, m, m} such that In— —m ■ fn — 1. Thus the 



covariant one-forms can be written in the Kinnersley-like forms |L5| as 
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I = dv — a sin 2 Odtp , 

n = ^ (dv — a sin 2 Odcpj — dr , 

m = -j=-^ {isin^ [adv — (r 2 + a 2 )dip] — Y,d9} , 

ra = -J^ {—i sin 9 [adv — (r 2 + a 2 )d(p] — Yid9} , 

and their corresponding directional derivatives are 

^ dr ' 



a _ r 2 +a 2 d _|_ A_d_ _|_ a_d_ 
— S dv "T" 2£ <9r "T" S 9^ ' 

5 =7fe( zasin ^ + i + 



_J §_ 

sin <9</3 



Inserting for the following relations among spin-coefficients |T4| 



A»-7 
7r — a 

0-T 



— r 
S ' 




r-M 


ia cos eA 


2£ 


2S 2 ' 


cote 


ira sin e 


2^,9 




cote 


a 2 sin e cos 


2 v / 2p* 





into four coupled Chandrasekhar equations [1 L| in the Newman- Penrose formalism [14 



(£ + e - ~ p )F x + (6 + ft- - a)F 2 = i$G x , 
( A + n - 7 )F 2 + (5 + (3 - t)F x = l -f 2 G 2 , 
(D + e* — p*)G 2 -{5 + tt*- a*)G x = ^§F 2 , 
(A + fi* - r)Gi + r*)G 2 = ^F x , 

where //q is the mass of Dirac particles, one can obtain 

Id |_ r \ p j_ 1 ( r ira si 

- [w + sj + 7T P \ L s 



IP-^J^+^^ r 



F 2 = 


72 Gl 






Gi = 


^F 2 


G 2 = 





here we have defined operators 
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r-M + 2a— + 2(r 2 + a 2 ) — 
dip dv 



or 

d 1 . id . . . d 

C = ^+o cot ^ — zasmtf— , 

dv 2 sin 9 dip dv 

d I id 

£ = ^ + -cot6»+^— — + iasin6»— . 
oa 2 sin 9 dip dv 



By substituting 



Pi, F 2 



P2 , G\ 



Qi j G2 



Q2, 



into Eq. (El), they have the form 



-I: Fx + £P 2 = inoffQi , ADP 2 + £t Pl = i^ p *Q 2 



(7) 



III. EVENT HORIZON 



Eq. (|7|) can not be decoupled except in the case of Kerr black hole fill] (M = const) 
or in the case of Vaidya space-time ]TB| (a = 0). However, to deal with the problem of 
Hawking radiation, one may concern about the behavior of Eq. (|7|) near the horizon only. 
As the space-time we consider at present is symmetric about the 99-axis, one can introduce 
a GTCT as does in Ref. |§ 

r * = r + ±\n{r-r H (v,6)], 

= v — v , 9* = 9 — 9 , 

where rn = rn{v, 9) is the location of the event horizon, and k is an adjustable parameter. 
All parameters k, Vq and #0 are constant under the tortoise transformation. The derivatives 
change correspondingly as 



d_ 

Or 


_ _d_ 

dr* 


+ 1 9 

2«() — m) dr* 


d_ 

dv 


- -JL 

dv* 


r H d 
2n(r— rn) dr* 


d 

de 


d 

de, 


r H d 
2re(r— rjj) dr* 



where f # = is the rate of the event horizon varying in time, r' H = -j^j- is its rate changing 
with the angle 9. 
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Under the transformations (§), Eq. (0) can be reduced to 
£-P x + (r' H - ia sin 8 f H ) £P 2 = , 

- (r' H + iasmd f H ) -§-P x + [A H - 2(r% + a 2 )r H ] -§-P 2 = , 



(9) 



and 



{r' H + ia sin Q r H ) £rQi ~ 4rQi = 



fr* ^ 1 dr* 

[A H - 2(r| + a 2 )r„] ^Qi + K - *a sin £ r H ) ^-Q 2 = 



(10) 



after being taken the r — > r H (v , 9 ), v — > v and 9 ^ 6 limits. We have denoted A H = 
r 2 H - 2M(v)r H + a 2 . 

If the derivatives -£rPi, ~£rP2, ~£~Qi an d ~£rQ2 in Eqs. (PJIHQ are nonzero, the existence 
condition of non-trial solutions for Pi, P 2 , Q± and Q 2 is that the determinant of Eqs. (|9| JT0|) 
vanishes, which gives the following equation to determine the location of horizon 

A H - 2(r 2 H + a 2 )r H + a 2 sin 2 6 f 2 H + r' H 2 = , (11) 

Eq. flTT|) just is the equation which can be inferred from the null-surface condition 
g^difdjf = 0, and f(v, r, 6) = 0, namely r = r(v, 6). 

The location of the event horizon is in accord with that obtained in the case of Klein- 
Gordon field equation ||. It means that the location of the horizon is shown as 

M±[M 2 - (a 2 sin 2 6 r 2 H + r' H 2 )(l - 2r H ) - a 2 (l - 2r H f} 1/2 n 

1 - 2r H 



IV. HAWKING TEMPERATURE 



Apparently one can find that the Chandrasekhar-Dirac equations (0) could be satisfied 
by identifying Q\, Q 2 with P 2 *, —P*, respectively. So he may deal with a pair of components 
Pi, P 2 only. Now, we rewrite the Dirac equation in a two-component form as 



.JL 

dr 



c 

£t AV 



\ P2 J 



\ Q2 J 



' -AV £HQi' 



V 



^ dr ) 



\ Q2 J 



-IfJLQp 



\ P2 J 



(13) 



and obtain its corresponding second-order form for the two-component spinor (Pi,P 2 ) as 
follows 
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AD J: + ££t + J, ( m sin #£t _ A £)] P 

+ { A (£ + z a S in 0£>) + za sin [2Mr f + M] } P 2 = /xgEPi , 
§-AV + £t£ + i (^a sin 0£ - AP)1 P 2 - i (£t + j a sin 0£) P x = /i 2 £P 2 , 



(14) 



Given the GTCT in Eq. (||) and after some tedious calculations, the limiting form of 
Eq. (|H]) , when r approaches r#(f , 6*0), f goes to v and 6* goes to 9 Q , reads 

- 2r H ) - M 



2 A H - 2r H (r 2 H + a 2 ) 
2 



+2(r 2 H + a 2 - r H a 2 sin 2 9 



02 o2 

—P 1+ 2a(l-r ff )— — Pi 



d 2 

p i " 2r H7r^7rA - {^(1 - 3r^) - M 



dr*dv* dr^dO* 

A H - r H (r 2 H + a 2 ) + ia sin 9 r' H d 

P* H dr * 1 



+ ^2iMr H a sin 9 - 
and 

> H (l-2r H )-M 



A h [t' h — ia sin 9 (r H + 1)] + 2i a sin 9 r H (r H + a 2 )) d 

P*H J 



A = 0, (15) 



K 



+ 2A H -2r H (r 2 H + a 2 ) 



+2a(l - r H ) 
4 



^2 ^2 

P> - 2r'„ - — P 2 - {M - r H (l + rjr) + r" H + f H a 2 sin 2 O + ^ cot O 



dr*dip~ * ' H dr*d9* 
Ah — rH{r 2 H + a 2 ) + ia sin 9qt' h d r' H + ia sin 9 (r H ~ I) d 



^ dr* ^ 2 



P = . 



(16) 



Ph ' ar * P*h d r * 

where p* H = rn + ia cos 9$. In the calculations, we have used the L'Hosptial's rule to treat 
an infinite form of jj-type. 

One can adjust the parameter k such that it satisfies 

r H (l - 2t h ) ~ M 



H 



+ 2A^ - 2r H (r H + a 2 ) = r 2 H + a 2 - r H a 2 sin 2 9 , 



(17) 



which means the temperature of the horizon is 

r H (l - 2r H ) - M 



" 2 H + a 2 - r H a 2 sin 2 9 - 2A H + 2r H (r 2 H + 



r H (l - 2r H ) - M 



" (r% + a 2 - r H a 2 sin 2 O )(1 - 2r H ) + 2r^ 2 ' 

where we have used Eq. (p"T|). 
Using Eq. (g), namely 

A Pl = _( r ' -i a8hi Q Q f H )JLp 2 8 P z = r H + ^^^r H d 
dr„ 1 H dr* 2 ' dr„ 2 A H - 2{r 2 H + a 2 ) dr* 



(18) 



Pi 
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Eq. (|15D and Eq. (|T^) can be recast into the standard wave equation near the horizon in an 
united form 



£)2 f)2 C)2 f)2 f) 

— M> + + 2{l——y + 2C 3 —-^ + 2(C 2 + id)— # = , 

or: or*ov* or*0(p or*o6* or* 



dr*dv* 



(19) 



where Q is the angular velocity of the event horizon of the evaporating Kerr black hole, 

a(l - r H ) 



Q = 



r 2 H + a 2 — f#a 2 sin 2 6 ' 



C, 



— r 



r 2 H + a 2 — r#a 2 sin 2 #o 



while both C\ and C*2 are real constants, 

-1 



Co 



2(r 2 H + a 2 - r H o? sin 2 



[r H (l-4r H ) -M + r' H 



... 2 • 2 q , / + /) 2Mr H f H a 2 sin 2 # 
+rjya sin # + r H cot £y + 



A H - 2r H {r 2 H + a 2 ] 



Ci 



2(r 2 H + a 2 - r H a 2 sin 2 9 ) 



-riia cos + 



2Mrnr ! H a sin #0 
A ff -2r fl (r| + a 2 ; 



for ^ = Pi, and 



for \l> = P 2 - 



C 2 
Ci 



M - r H + r" H + r^a 2 sin 2 6> + cot C 
2(r|^ + a 2 — r#a 2 sin 2 # ) 
ri/a cos 9 
2{r 2 H + a 2 — r H a 2 sin 2 9 ) 



V. THERMAL RADIATION SPECTRUM 



Now separating variables as ^ = R(r*)®(Q*)e %<ymip wv *> , one has 

6' = A6 , 

R" + 2(C + %C X + imtt - iuj)R' = , 

where A is a real constant introduced in the separation of variables, Cq 
solutions are 

= e Ae * , 

R rv, e 2i(w-mn-C*i)r*-2C*or, . _ 



(20) 



XC 3 + C 2 . The 



(21) 
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The ingoing wave solution and the outgoing wave solution to Eq. (|19"|) are, respectively, 



(22) 

The outgoing wave \l/out is not analytic at the event horizon r = r#, but can be analyti- 
cally continued from the outside of the hole into the inside of the hole by the lower complex 
r-plane, 

(r - r H ) -> (r# - r)e" i7r 

to 

{fj ^ __ e -iaj^»+im(/3+Ae, e 2i(aj-mf2-Ci)r,-2Cor« e i7rCo/Kg7r(a;-mSl-Ci)/K ^ (23) 

The relative scattering probability at the event horizon is 

g~27r(u;— mil- Ci)/k (24) 



^out 

Following the method suggested by Damour and Ruffini 0], and extended by Sannan ||, 
the Fermionic spectrum of Hawking radiation of Dirac particles from the black hole is easily 
obtained 



VI. CONCLUSION 

Equations (|l"2"D and (|1|) give the location and the temperature of event horizon of the 
variable-mass Kerr black hole, which depend not only on the advanced time v but also on the 
angle 8. They can recover the well-known results previously obtained by others. Eq. (|25|) 
shows the thermal radiation spectrum of Dirac particles in the non-stationary Kerr space- 
time, in which a residual term C\ appears. The difference between Bosonic spectrum and 
Fermionic spectrum is that C\ is absent in the spectrum of Klein-Gordon scalar particles. 
Also C\ vanishes when a black hole is stationary (fjj — r' H — 0) or it has a zero angular 
momentum (a = 0). Thus this new effect maybe arise from the interaction between the spin 
of Dirac particle and the evaporating black hole. 

In conclusion, we succeed in dealing with Chandrasekhar-Dirac equation in the non- 
stationary Kerr black hole. Under the generalized tortoise coordinate transformation, each 
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second-order equation induced from of Chandrasekhar-Dirac equation takes the standard 
form of wave equation near the event horizon, to which separation of variables is possible. A 
key step in our arguments is that we have considered the asymptotic form of the first-order 
equation near the event horizon, and obtained formulas relating the first-order derivatives 
of different component which enable us to eliminate the crossing-terms in the second-order 
equation. The location and the temperature of event horizon of the variable-mass Kerr 
black hole are just the same as that obtained in the discussing on thermal radiation of 
Klein-Gordon particles in the same space-time. But the spectrum of Dirac particles displays 
another new effect which is absent from the spectrum of Klein-Gordon particles. This effect 
is a kind of spin-rotation coupling effect. 



APPENDIX: NEWMAN-PENROSE COEFFICIENTS 

In this appendix, we make use of the method of exterior differentiation to calculate 
the Newman-Penrose coefficients in the non-stationary Kerr space-time. The complex null- 
tetrad {I, n, m, Tn} satisfies the orthogonal conditions In— —m • rn — 1. We choose the 
covariant one-forms similar to the Kinnersley forms as 



(Al) 



I = dv — a sin 2 0d<p , 
n — A (dv — a sin 2 0d<pj — dr , 
m = -j=-^ {i sin 9 [adv — (r 2 + a 2 )d^} — Y,d6} , 
rn =-j=-{-i sin [adv - (r 2 + a 2 )dip] - Y,d6} . 

The exterior differentiation of the one-forms bases gives 
dl = v^smecose (m + np) A i-?iaime rnA rn, 

J„ / rA r—M \ a i iMra sin 9 cos ( TTl TTl \ a / ia cos 9 A a -^r 

an - ^ S 2 2S J' iAl 2V2S \p p* ) £ 2 mAm > 

dm = ^0^n A I - i(n - £l) Am + ^fcote + ^jmAra, 
dm = - ^y' Tt A I - ± (n - A;) a m - ^ (cot 6 - m Am . 

In general, the common relations among the exterior differentiations of one-forms bases 
hold: 



(A2) 
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dl = (e + e*)n A Z + (a + /3* - r*)Z A m + (a* + (3 - r)i A m 

— R*n Am — kn A m + (p — p*)m A m , 

c?n = (7 + 7*)n A Z + z/Z A m + z/*Z A m 

+ (7r — a — /3*)n A m + (n* — a* — B)n A m + — u*)m A m , 

_ (A3) 

dm = (ff* + r)n A Z + (7 - 7* + /x*)Z A m + A*Z A m 

+ (e — e* — p)n Am — an A m + (/9 — a*)m A m , 
c?m = (-7T + r*)n A I + AZ A m + (7* - 7 + /z)Z A m 

- cr*n A m + (e* — e — /5*)n A m + (a - /3*)m A m . 

where a star * stands for complex conjugate as usual. 



It is not difficult to determine the twelve Newman-Penrose complex coefficients |14 
the above null-tetrad as follows 

0, 

rA r-M 



111 



K = 


0. 


A = 


0, 


0" 


P = 


1 


e = 


ia cos 9 


7 


P * ' 


S 


P = 


A 

2Sp* ' 


T = 


iasin(9 


;y 


a = 


7T-/?*, 


7T = 


ia sin 


, /3 





2E2 2S ' (A4) 

Mria sin 
cot _|_ ira sin 9 



2 v / 2p* v^Sp* ' 
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